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ABSTRACT 

We present projected rotational velocity values for 97 Galactic, 55 SMC, and 
106 LMC 0-B type stars from archival FUSE observations. The evolved and un- 
evolved samples from each environment are compared through the Kolmogorov- 
Smirnov test to determine if the distribution of equatorial rotational velocities is 
metallicity dependent for these massive objects. Stellar interior models predict 
that massive stars with SMC metallicity will have significantly reduced angular 
momentum loss on the main sequence compared to their Galactic counterparts. 
Our results find some support for this prediction but also show that even at 
Galactic metallicity, evolved and unevolved massive stars have fairly similar frac- 
tions of stars with large V sin i values. Macroturbulent broadening that is present 
in the spectral features of Galactic evolved massive stars is lower in the LMC and 
SMC samples. This suggests the processes that lead to macroturbulence are de- 
pendent upon metalhcity. 

Subject headings: Stars: spectroscopic — stars: early-type — stars: rotation — 
Magellanic Clouds — ultraviolet: stars 

1. Introduction 



Stellar interior models are the backbone of modern stellar astrophysics. Their predic- 
tions are used in a myriad of applications including determinations of stellar masses from 
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observed luminosities and temperatures, studies of star clusters, population synthesis, stellar 
nucleosynthesis, chemical evolution, etc. Many factors are included in these extremely com- 
plex interior codes: overshooting, mass loss rates, opacities, convection. The latest addition 
to t hese models has been the inclusion of rotation and its subsequent effects on interior mix- 



ing (IMeynet fc Maeder I l2000l : iHeger fc Langer 1 120001 : ?). With each additional advance in 



these codes, observers have been tasked with testing the new models by finding observable 
diagnostics. What do the new models predict that can be observed to either agree or disagree 
with the model? For massive stars the newest models are very interesting. Rotation induces 
interior mixing that changes the luminosities, lifetimes, effective temperatures, and surface 
abundances of these massive stars. The combined effect of mass loss and angular momentum 
in the models results in angular momentum loss; as the mass comes off it takes with it any 
angular momentum it may have had. This results in an evolution of surface rotation speed 
that varies depending upon the amount of mass loss. 

The most striking prediction of the new models is the effect of rotation on the internal 
structure and subsequent luminosity and main sequence lifetime of the star. The models 
with rotation predict evolutionary tracks that become more luminous. This results from two 
effects: rotational mixing of H-fuel into the convective core and transport of He and other 
H-burning by-products into the radiative envelope. The first effect incr eases the mass of 



the co nvective core, while the second lowers the opacity of the envelope. iHeger &: Langer 



(I2OOOI ) state, "As a result of these dependencies of the stellar evolutionary tracks in the HR 
diagram on the initial rotation rate, a given point in the HR diagram is not uniquely related 
to a single initial mass, even for core hydrogen burning stars." This introduces a serious 
scatter in the mass-luminosity relation. In order to determine the evolutionary mass for an 
object, the luminosity, effective temperature, and the initial equatorial velocity of the star are 
needed, but this last parameter is indeterminable. One can attempt to test this prediction 
of rotationally induced overluminosity. To test if a single star is overluminous for its mass, 
first requires knowledge of the mass of the star. The only direct means to determine the 
mass is from binary systems. There are cases where members of a binary system with large 



equat orial velocities do appear overluminous for their masses (jPenny. Gies. &: Bagnuolo 



19991 ). However other effects can also cause binary stars to appear overluminous, mass 
transfer being the most common. Also members of a binary pair exchange spin and orbital 
angular momentum as they evolve due to tidal interactions, making them poor tests of single 
star models. 

The models also predict how the surface rotation rates should change as the stars evolve. 
For solar metaUicity stars {Zmw = 0.020) the most massive should slow their rotation 
drastically while still on the main sequence (MS). There are two competing factors that 
govern the surface rotation rates: (1) transfer of angular momentum from the inner parts to 
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the outer ones through meridional circulation and (2) angular momentum loss at the surface 
of the star through mass loss. The first of these is present in all rotating stars and is primarily 
dependent upon the rotation rate: the faster the spin, the more mixing occurs. The mass loss 
is from large stellar winds that are driven by line opacities associated with metallic elements. 
For stars at lower metallicities (i.e., Small Magellanic Cloud, Zsmc = 0.004 and Large 
Magellanic Cloud, Zlmc = 0.007), there are significantly lower numbers of metal atoms and 
therefore lower line opacity. Mass loss rates are lower and subsequent angular momentum 
loss rates are lower too. The rotation rates of massive stars at lower Z are predicted to 
remain almost constant during the MS. For example Figure 10 of iMeynet &: Maeder I (120001 ) 
shows the evolution of surface rotational velocities as a function of time for a 20 Mq star at 
solar and SMC metallicity with an initial rotation speed of 300 km s~^. By the late main 
sequence the rotation rate of the Galactic star has decreased to 140 km s~^, while its SMC 
counterpart still has a rotation rate of 220 km s~^. 

It is interesting to consider what spectral class stars these models refer to, as stars do 
not come labeled with their initial masses. The new models predict that stars with initial 
masses smaller than 15 Mq do not slow significantly while on the MS even in the Milky 
Way. There is no metallicity dependence on their rotational evolution, i.e., a 12 Mq star 
at both solar and SMC metallicity have similar behavior. Our prime targets for testing 
purposes should be stars that started with 15 or more solar masses. In addition we see that 
after the end of the MS, the surface rotational velocities of even these massive stars drop 
precipitously in all metallicity environments. A star with an initial mass larger than 15 Mq 
spends its MS life as an 0-type star (see Fig. 1). In fact, we see that despite expectations 
that luminosity classes of giant or supergiant describe post-MS stars, this is not possible for 
O spectral types. All 0-type stars, regardless of their luminosity class, are core-hydrogen 
burning objects. They will eventually become B and A supergiants, but when they do they 
will have greatly increased radii, and decreased surface rotational velocities. They also do 
not exclusively populate these types. In Figure 1 we over plot the effective tera perature line 
of 14900K, which corresponds to the luminosity class B3 I (lBohm-Viteiiselll98lf) . Even a star 

Teff = 15100K, 



1992 



with initial mass of 7 Mq becomes a B3 I (according to ISchaller et al. 
\ogg = 3.25 occurs at age =4.32 x 10'' yrs for the 7 Mq track). Early B-type supergiants 
are not an exclusive group and for stars with Mj„j > I5M0, they represent stars past the 
end of the MS. As such they are not optimal targets for testing the predictions from the new 
stellar interior models. 

There exists observational evidence that the predictions for solar metallicity stars are 
correct. Galactic 0-type stars do slow down as they evolve along the main sequence 
( jConti fc EbbetsI 119771 : iPennyl Il996l : iHowarth et al.l 119971 ) . There are few broad-lined stars 
among the more massive and more evolved objects. The fastest rotators are found among 
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the lower mass O stars. This agrees well with predictions: that the rotation rates should 
decrease while still on the MS; that the larger the initial mass, the faster the decrease in 
rotation; and that the more evolved MS stars should show the biggest decrease in rotation. 
Observationally only V sini (projected rotational velocity) is determined in these studies due 
to the orientation of the polar axis to the line of sight (i). But with a large enough sample 
we can assume a random distribution of this inclination angle. While it is possible that some 
massive and evolved objects are observed pole-on, with a large sample it is unlikely that 
orientation effects can explain the observed trends. 

Indirect support that massi ve stars at low metallcity do NOT slow down comes from 
a study by iMaeder et al.l (119991 ) who found a significant increase in the number ratio of 
Be/(B-|-Be) for decreasing metallicities in Galactic and MC clusters. They argue that the 
Be phenomenon is closely linked to rapid rotation, and so if the massive stars with low Z 
do not slow down, perhaps they produce more Be stars. We should note here that a Be star 
is a non-supergiant B-type star that has a significantly lower initial mass than the 0-type 
objects. The models produced for low Z massive stars do not predict that stars with initial 
masses relating to Be-type will have metallicity dependent decrease in rotation rates while 
on the MS. However the initial velocities of the se stars may vary w i th Z. The initial rotation 
rates of massive stars at low Z are not known. iMaeder fc Meynetl (120011 ) have assumed that 
the initial rates are the same in the MC and Galaxy (since observational data are lacking). 
But the star formation process itself may be altered at smaller metallicity. Although no 
numerical models of the formation of massive stars at low Z have been produced, there are 
some possible origins of higher initial rotational velocities. A lower Z environment would 
have lower dust content and fewer metallic ions present in star forming regions so that the 
magnetic field of the contracting central mass is less coupled to the surrounding region. Also 
the lifetime of the accretion disk is likely related to its opacity, which will decrease with lower 
metallicity. These factors would result in less angular momentum loss during star formation, 
leading to higher initial rotational velocities. 

There are few direct measurements of the rotation rates of massive stars at low metallic- 



ities. IMaeder et al.l (119991 ) state "Of course, direct observations of V sini in LMC and SMC 
clusters are very much needed in order to substantiate the above results." Only a few studies 
have been published to date on the rotation rates of massive stars a t low metallicity. Some 
have focused on early B type stars (lKellerll200ll : iDufton et al.ll2006l ). but again we caution 
that these cooler objects are not exclusively produced by stars with 15 or more solar masses 
and are often significantly larger than their terminal age main sequence (TAMS) radii. 

We mad e a study of 44 -type stars in the SMC and LMC utilizing archival HST 
observations (jPenny et al.l 120041 ). We compared the cumulative distribution function (CDF) 
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of the V sin i values for Galactic, LMC and SMC evolved (giant, bright giant, and supergiant) 
stars. The CDF gives the fraction of stars having a V sin i less than a specific upper value 
that ranges from zero to the maximum Vsini observed. The surprising result is that the 
distributions are very similar, with a slightly larger fraction of LMC stars rotating faster, 
but a slightly smaller fraction of SMC stars, with lowest Z, rotating slower than the Galactic 
stars. The two sample Kolmolgorov-Smirnov (K-S) test compares two independent samples 
and determines a statistic, p describing the similarity of the CDFs. For low values of p {< 
5 — 10%), the null hypothesis that the two samples are drawn from the same parent population 
is rejected. Comparing the LMC and SMC samples with Galactic stars, their K-S p are 44% 
and 42%, respectively, indicating that they are drawn from the same parent population as 
the Galactic sample. Thus, there is no visible trend of CDFs with metallicity. The sample of 
unevolved (dwarf and subgiant) stars was too small to make a reliable comparison, but the 
average projected rotational velocities in each group appear to be comparable: < V sini >= 
131 ± 93, 117 ± 31, and 78 ± 45 km s~^ for the Galactic, LMC, and SMC unevolved stars, 
respectively. The results of this study should be treated cautiously for two reasons. First, 
the samples of LMC and SMC targets are not large enough to assume a completely random 
distribution of inclination angle on the sky. Second, 18 of the HST/STIS SMC observations 
were obtained through one observing program. This program's goal was to study the UV 
characteristics of massive stars at low Z. These targets may have a bias toward small projected 
rotational velocities, in order to better study the line depths as a function of metallicity. 



Hillier et al.l (120031 ) and iBouret et al.l (120031 ) examined a total of 8 0-type stars in the 



SMC by fi tting rotationally bro adened synthetic spectra from the atmosphere model code 
TLUSTY ( iHubeny fc Lang|l995l ) to the observed spectra. In all but one case the determined 
projected rotational velocities were < 100 km s~^, much smaller than the predicted equatorial 
rotation speeds of ~ 200 km s~^. There are several possible explanations: (1) the initial 
rotation rates may not be as large as 300 km s~^ for massive stars in the SMC; (2) the stars 
may experience more spin down than predicted; and (3) with only 8 targets, the sample is 
just too small considering the random distribution of inclination. 

The most recent survey of rotational velocities of 0-type stars at low Z is a study by 
( 20061) . which is part of the larger VLT-FLAMES survey of massive stars 
20051 ). They determined Vsini values for 28 0-type stars in the SMC. A 



Mokiem et a. 



(lEvans et al 



comparison of their unevolved targets (classes V, IV) to evolved targets (classes I, II, & 
III) shows that the unevolved sample contains relatively more fast rotators tha n the evolved 



stars. They point o ut that this trend is also seen in Galactic O type stars (jPennyl Il996 



Howarth et al.l 119971 ). New stellar interior models predict this behavior for solar metallicity 
stars but not for massive stars at Zsmc = 0.00 4. They also compare their 21 unevolved 
SMC stars to 66 unevolved Galactic stars from iPennyl (119961 ). The K-S p is only 13%, 
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suggesting that there is a possible difference i n the initial rotation rates of massive stars due 



to metallicity. However, iMokiem et al.l (120061 ) warn that because of their small sample size, 



this result should be verified with a larger SMC sample. 

All of the previous V^sini surveys of massive stars at low Z suffer from one common 
major fault, a small sample size. Here we present the results of our project to test the 
treatment of angular momentum in the new stellar interior models through a large scale 
survey of projected rotational velocities of 0-type stars in three metallicity environments: 
the Milky Way {Zmw = 0.020), LMC {Zlmc = 0.007), and SMC (Zsmc = 0.004). The Far 
Ultraviolet Spectrographic Explorer (FUSE) archive at the Multimission Archive at Space 
Telescope (MAST) contains spectra of 161 LMC and SMC stars with spectral classes between 
B2 - 02. The targets are 120 evolved (luminosity classes I, II, & HI) and 41 unevolved (IV 
& V luminosity classes) stars in these low Z environments. In addition MAST contains 
spectra of approximately 100 Galactic 0-type stars with known Vsini values. This is a 
large sample of 0-type spectra, all observed with the same instrument. It contains both 
cluster and field stars. The observations were not part of any single campaign, in fact they 
were obtained through 38 different observing programs. These stars represent a completely 
uniform, unbiased sample. Comparisons of the Vsini distributions from the unevolved 
MC stars to similar classed Galactic stars will determine whether the initial rotation rates 
of massive stars have a metallicity dependence. Comparisons of the evolved stars in each 
environment to their unevolved counterparts will show whether the rotation rates are age 
dependent at that metallicity. 



2. Projected Rotational Velocities from UV lines 

For the extremely massive 0-type st ars, the major na e chani sms that affect stellar line 



widths are rotation and macroturbulence. IConti fc Ebbetd (119771 ) completed the first large- 



scale survey of rotational velocities of Galactic 0-type stars. They noted in their sample 
that no 0-type star displayed extremely sharp lines that would be expected with a pole-on 
orientation and suggested that large-scale atmospheric motions contribute to the broadening 
of the spectral features. In determining their V sin i estimates they adopted the convention 
that only rotational broadening was present. They were awar e this is not in general true 



and note, as do others who have adopted a similar convention (jPennylll996l : iHowarth et al. 



19971 ). that their quoted Vsini is actually a parameter which is closely related to, but not 



necessarily, the exact projected equatorial velocities of the stars. Authors who have investi- 



(Rvans et al. 


2002; 


Dufton et al. 


2006) 
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appea r to be dependent u pon metallicity as both the Galactic (IRyans et aLll2002l ) and the 
SMC ( iDufton et al.l l2006l ) samples result in similar values for macroturbulent broadening. 
While macroturbulence is undoubtedly a factor in the line widths of 0-type stars, we can 
make a statistical argument that if each sample (Galactic, SMC, LMC) has the same contri- 
bution to line widths from macroturbulence, then the distribution of Vsini from each can 
be compared for differences in equatorial rotational velocities (see discussion in §5). 



Our previous studies of projected rotational velocities (|Pennylll996l : iPenny et al.ll2004l ) 
utilized archival high-resolution spectra from the International Ultraviolet Explorer (lUE) 
and the Space Telescope Imaging Spectrograph (STIS) in the wavelength region 1200 - 
1900 A. The ultraviolet is actually a preferable region for determining projected rotational 
velocities for very massive stars. Unlike the photospheric lines seen in optical spectra, the 
UV lines are formed deeper in the star and are usually less contaminated by wind effects. 
Also there are thousands of lines present in the UV, compared to tens in the optical. The 
methodology used to obtain V sin i values from UV spectra differs from that used for optical 
spectra. We cross-correlate the spectrum of a narrow-lined star with that of a test star, 
which results in a cross-correlation function (CCF) that represents a "superline" of the test 
star. The superline will have an observed width related to the line width of the test star and 
that of the narrow-lined star. In the case of a double-lined binary, the CCF may display two 
peaks, which are separated by the radial velocity difference of the two stars. A calibration 
between CCF width and V sin i is obtained by cross-correlating a sample of stars with known 
projected rotational velocities. 



3. Observations and Reductions 

The CCF method can also be utilized with FUSE observations of 0-type stars. We 
obtained FUSE archival spectra of 97 Galactic, 55 SMC, and 106 LMC massive stars that 
are available from MAST. For each star we obtain the calibrated (CALFUSE) exposure 
files from the LiF2 channel and detector segment A (hereafter LiF2A spectra). The LiF2A 
spectra cover only the wavelength region of 1086 - 1183 A. Below this region the spectra 
are severely affected by interstellar features and are not suitable for our purpose. Each ob- 
servation is actually composed of several exposures. The exposure spectra were extracted 
using the CALFUSE routine MRDFITS. Following this basic reduction we use a series of 
routines written in the Interactive Data Languagqj to further reduce the spectra. The ex- 
tracted exposure spectra then were coadded with appropriate wavelength shifts, determined 



^IDL is a registered trademark of ITT Corporation 
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by cross-correlating the exposure with the highest flux value (reference spectrum) to all 
others, creating a weighted average spectrum based upon exposure lengths for each obser- 
vation. These steps are consolidated into one routine COMBINE that writes the coadded 
spectrum to a file. Data spikes are removed in SPIKEOUT, which replaces any bad pixels 
with a median of surrounding pixels. One of the most important steps in this reduction 
is the placement of the data on a log A grid, in which each pixel step corresponds to a 
uniform velocity step of 1 pixel = 10 km s~^. This is done in the routine BINFUSE. The 
spectroscopic resolving power of the raw FUSE spectra is 20,000 ± 2000, however we further 
degrade this value by smoothing the data by a Gaussian transfer function (FWHM = 40 
km s~^) in GSMOOTH. Another important step is the alignment of all spectra of a star 
on the same wavelength system. In order to assure that all spectra of a given star are not 
shifted with respect to each other, the interstellar lines present (assumed to have a constant 
velocity) are used to align the system spectra. ISMALIGNFUSE does this by creating an 
average spectrum, cross-correlating the regions containing strong interstellar lines (which are 
picked interactively), and shifting each spectrum to an average wavelength grid. ISMALIGN- 
FUSE also removes not only the strong interstellar features, but also those chosen (again, 
interactively) to be weak interstellar features. We cannot overemphasize the importance of 
removing the cons iderable number o f inte rs tellar lines, whose id entific ation was greatly aide d 
by earlier work of iBarnstedt et all kood \. IPellerin et aD J2002h . and IWalborn et al.l J2002[ ). 
The spectra are then rectified (routine FUSERECT) to a unit continuum, by fitting a spline 
function to a series of pseudo-continuum zones in an average spectrum and then dividing 
each of the spectra by the fit. The final spectral stacks are written to a file in WTSTKFUSE. 
All of these procedures are incorporated in a calling routine named MKSTKFUSE. 



4. Projected Rotational Velocities from Cross-Correlation Functions 

Target star spectra are cross-correlated with that of a relatively narrow-lined template 
star. Cross- correlation was done only in the wavelength regions 1100 - 1108 A and 1130 - 1165 
A (1139 - 1165 A, for stars with emission in the P V/ Si IV A1128 line), a region relatively 
free of wind features. While most of the lines in these regions are too weak and blended to 
identify individually, a few of the stronger features are O III AA1149.6, 1150.9, 1153.8 A, and 
S IV Al 138.2 A. The spectra are padded with 1000 km s~^ of artificial continuum on both 
ends to include the entire observed range in the cross-correlation. The CCF is the sum of 
the square of the differences between the test spectrum and the reference spectrum shifted 
in velocity from —1000 km s~^ to +1000 km s~^ at 10 km s~^ intervals. The functions are 
then rectified, inverted, and Gaussian fitted to obtain an estimate of the full-width at half- 
maximum (FWHM). In Figure 2 we present fitted CCFs for 4 Galactic targets, HD 46202, 
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HD 152233, HD 96917, and HD 41161 that were all made with the template spectrum of 
star AV 327 (09.7 Il-Ibw). We note that the 'bumps' present in the CCF of HD 41161 are 
typical for massive stars with large projected rotational velocities and are possible indicators 
of non-radial pulsations present in the star. 

The next step is developing a calibration of CCF width to V sin i usin g the Galactic 
targets. These stars have known Vsini values from iHowarth et al.l (119971 ). We created 
calibrations for two template stars AV 327 (09.5 Il-Ibw) and SK -66 100 (06 11(f)) in order 
to cover the large range in spectral class of our targets. These are presented in Figure 3. 
Because of the early type of SK -66 100, seven of the latest type stars did not produce 
significant CCF peaks and were not used to create the calibration for this template. In 
Table 1, we lis t the Galactic targets along with their spectral classes, Vsini values from 
Howarth et al.l (119971 ). Vsini values computed from the calibrati ons of each ten i plate star, 
and the FUSE observations used. Spectral classes are as listed in iHowarth et al.l (119971 ). 



For each target from the Magellanic Clouds, we follow a similar procedure. All targets 
are cross-correlated with both templates to determine which template creates the CCF with 
the highest peak. For target stars with more than one spectrum, each spectrum's CCF is 
examined to determine if the peak is stationary. A CCF peak that moves from observation 
to observation is indicative of a single-lined spectroscopic binary, as in the case of NGC346- 
MPG324 which shows a peak shift of 110 km s~^ between two observations. In these cases, 
each CCF is fit individually. For targets with more than one observation and stationary 
peaks, the average CCF is fitted. A few of the stars' CCFs had double peaks, indicating that 
they were double-lined spectroscopic binaries. These CCFs were fit with a double Gaussian. 
Then using the calibrations developed above a V sin i value is determined from the Gaussian 
width of each fitted CCF. In Tables 2 and 3 we list the SMC and LMC targets, their spectral 
classifications, V sin i values from this paper, along with any previous values, the template 
star utilized, FUSE observations, and an additional column for other comments. For those 
CCF with two peaks the weaker pe ak's V sin i is g i ven in pare nthes es. In Figure 4 we plot 
our V sin i values versus those from iMokiem et al.l (l2006l . 120071 ) and iPenny et al.l (120041 ) for 
those targets in common. From the scatter in Figure 4, we estimate that the typical errors 
for our ^sini values are 10 km s~^. 



5. Results and Conclusions 

We can compare the projected rotational velocity distributions of the different groups 
of stars by calculating their cumulative distribution function (CDF). Again, the CDF gives 
the fraction of stars having a V sin i less than a specific upper value that ranges from zero to 
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the maximum V sin i observed. Not all of our targets stars are used in compiling the CDF. 
Stars with a 'NO' in the final column of Tables 2 and 3 are excluded from our analysis for 
various reasons. For the B-type stars in our sample, we perform an additional analysis to 
determine: if they formed from stars with initial masses larger than 15 Mq, and the fraction 
of their current radii to their TAMS radii. For all the stars we also wish to insure that 
extremely close binary systems are not present. We calcu l ate M hni from their V, B — V 
values, the {B — V)o cali bration of IShuU &: van Steenbera (119851 ) . Ay = 3.1 x E{B — V), 
bolometric corrections of iHowarth fc Prinjal (119891 ). and a distance modulus of 18.3 and 
19.1, for the LMC and SMC re s pectiv ely. We adopt t he effective teraperatu res from the 
calibration of iHowarth fc Prinial (1989) for 0- t ypes and iBohm-Vitensd ( 119811 ) for B-types. 
From the evolutionary tracks of ISchaerer et al.l (Il993l ). we extrapolate initial masses, current 
masses, log g, current radii, and TAMS radii. As we expected a large number of our B-type 
targets are inappropriate for our analysis. None of the B-stars have initial masses smaller 
than 15M0. However, thirty-six targets (11 SMC and 25 LMC) are a factor of 50% or more 
larger than their TAMS radii. The lack of lower mass B-type supergiants almost certainly 
refiects their lower luminosities and the subsequent difficulty of obtaining FUSE spectra of 
those stars. For the entire FUSE sample we also compare the interpolated log g values 
with the luminosity class for discrepancies which might indicate that more than one star is 
present in the spectrum. F or o ur analysis w e also include 27 SMC and 26 LMC stars from 



Mokiem et al.l (l2006l . 120071 ) and lPenny et al.l (120041 ). We have already shown in Figure 4 that 



there are no systematic differences between the Vsini values from our FUSE analysis and 
those presented in these prior studies. These ad ditional stars a r e presented in Tables 4 and 



5, alo ng with the spectral classification given in iMokiem et al.l (l2006l . 120071 ) or iPenny et al. 
(I2OO4J ). and those authors' V sini values. 

We divide the data into the following samples: SMC unevolved (32 stars), SMC evolved 
(19 stars), LMC unevolved (36 stars), and LMC evolved (42 stars) . We also crea t e sara - 



ples for Galactic unevolved (79 stars) and evolved (56 s tars) from 



Howarth et al 



(11997t). 



Following the convention of ear lier studies ( Penny 1996; Howarth et al. 19971 : Penny et al. 



2OO4J : iMokiem et al.ll2006l . 120071 ). unevolved refers to stars with luminosity classes V & IV, 



while luminosity classes I & II are termed evolved. Class HI stars are omitted from both 
groups. Three stars are without explicit luminosity classes, AV 80, SK -67 266, and SK 
-69 279, and we assign these to the evolved category based on the notes from their listed 
spectral classification source. For each sample we create a CDF and plot these in compar- 
ison with other samples to determine the similarity. For each comparison a K-S statistic, 
D, and its corresponding probability, p, are calculated and these are presented in Table 6. 
The test statistic, D, is the maximum vertical variance between the two CDFs at any V sini 
value. The p value tells us the probability of obtaining a D value at least as extreme as 
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was observed, assuming that the two samples are drawn from the same parent population. 
So in cases where the p statistic is small, i.e., < 5%, we can say that the probability of 
finding such a large D between two samples that are drawn from identical populations is 
very small. Plots of the comparisons made are presented in Figures 5-7 and 9-10. At the 
suggestion of the referee, we determined the mass distribution of each of our samples. For 
the MC stars, we adopted the current evolutionary masses extrapolated from the Mm and 
Teflf analysis described in the paragra ph above. Mass estimates for the Galactic sample can 



be found in lHowarth fc Prinjal (119891 ). The mean masses of the samples are 37.8, 45.7, and 
34.4M0 for the unevolved stars of the Galaxy, LMC, and SMC, respectively, and the means 
are systematically higher among the evolved stars of the same three environments, 45.4, 47.4, 
and 4O.3M0, respectively. This difference is expected since the evolved 0-star samples will 
be biased towards higher masses (the advanced evolutionary stages of lower mass stars occur 
at cooler temperatures). The K-S test indicates that the three evolved sets have statistically 
indistinguishable mass distributions. On the other hand, the K-S test suggests that some 
differences may exist among the unevolved star mass distributions (particularly between the 
SMC and LMC samples where p = 0.7%), but given the relative small sample sizes, we doubt 
that these are significant. 

First we examine the evolved and unevolved Galactic samples to determine the level of 
difference we would expect in an environment where the stars do slow down as they evolve 
(Fig. 5). We see that the maximum variance in their CDF is 0.30 and this occurs at Vsini 
of 83 km s~^. This large a D value results in a p of 0.5%, but we should not be misled by this 
low value. At these smaller velocities the source of the difference between the evolved and 
unevolved samples is not a result of angular momentum loss, but from the larger amount 
of macroturbulence that is present in the evolved stars' photospheres, broadening their line 
profiles. At the higher velocities where we expect to see evidence of spin down, the largest 
divergence is ~ 0.12, which would give us a significantly higher p of 71%. This is not say 
that the evolved stars in the Galactic sample have not slowed down, but that this effect may 
be more subtle than we originally expected. 

In the LMC, the same comparison of unevolved to evolved stars has a very different 
appearance (Fig. 6). The CDF of both samples appear very similar, and the derived D = 0.15 
and corresponding p = 75% supports the null hypothesis that both are d rawn from the sam e 



parent distributi on. This resu lt is in agreement with the suggestion of IWolff et al.l (120081 ) 



using data from iHunter et al.l (120081 ). that massive stars in the LMC have similar Vsini 
distributions regardless of their logg values. Examining Figure 6, we notice that there is 
almost no divergence at low V sin i values, and a smaller variance at larger V sin i than we 
see in the Milky Way samples. The very high p value here is primarily due to the presence 
of slow rotators among both the unevolved and evolved samples, unlike the Galactic evolved 
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group. At high ^sini values, the CDFs of evolved to unevolved differ by ~ 0.09 which 
is only slightly smaller than the value for the Galactic samples. It is at the low velocity 
end of the CDF where the LMC samples differ from those of the Galactic samples. The 
good agreement between the low ends of the two CDFs for the LMC stars suggests that 
development of macroturbulence with evolution is not as large a factor in the photospheres 
of the massive stars in that environment and that the processes that lead to macroturbulent 
broadening may have a metallicity dependence. 

The CDFs of the same populations in the SMC are presented in Figure 7. Although 
the p value from the KS test is 23%, which would indicate that the samples are drawn from 
the same parent population, the maximum difference D is 0.29, which is very close to that 
from the Galactic samples. The larger p value from the similar D reflects the much smaller 
sample sizes in the SMC. At the high velocity end of the CDF, the fraction of evolved stars 
with V sini above 200 km s~^ is 11% compared to 13% for the unevolved. This does support 
the hypothesis that stars at SMC metallicity will not slow down as they evolve on the MS. 
At the other end of the CDF, there is a slightly larger fraction of unevolved stars with 
Vsini below f» 70 km s~^, but this difference is much smaller than we see in the Galactic 
comparison. The real divergence between the evolved and unevolved CDFs comes at the 
intervening velocities. The maximum divergence, 0.29, occurs at Vsini = 107 km s~^. The 
fractions of stars with V sin i at or below this value for t he evolved and unevol ved samples are 



0.84 and 0.55, respectively. This trend was also seen by lMokiem et al.l (120061 ). who surmised 
that the initial rotational velocity distribution in the SMC might vary significantly from that 
in the Galaxy. Why we see a disparity in the CDFs of the evolved and unevolved samples in 
between 80 — 190 km s~^ only in the SMC samples is not clear. We emphasize that the K-S 
test results accept the null hypothesis that the unevolved and evolved samples in the SMC 
are drawn from the same parent distribution. 

A primary purpose of this project is to examine observationally the amount of angular 
momentum loss during the MS lifetimes of massive stars and the effects of metallicity on 
this loss. Comparing the evolved to unevolved samples in Figures 5-7, we see that the D 
statistic at high Vsini values range from 0.12, 0.09, 0.02 in the Galaxy, LMC, and SMC, 
respectively. Taken by itself this is suggestive of a trend with decreasing Z. However, for the 
LMC and SMC comparisons, D statistics of 0.09 and 0.02 result in p values of 99.6% and 
100%, indicating there is no statistical difference between the evolved and unevolved CDFs 
at high Vsini. But this is also true in the Galaxy where we do expect to see the significant 
spin down between the evolved and unevolved samples. Here the corresponding p for a of 
0.12 is 71%, far above the 5% cutoff. Statistically the loss of angular momentum for these 
H-burning stars is not so different between these three metallicity environments. 
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Although our results may appear to cast doubt on the rotational evolution predicted 
by models, we suspect that the rather small changes we find are probably consistent with 
models. Unfortunately, at the moment there are no large grids available for rotational evolu- 
tio n that might be used to compare with our statistical results. The best set was calculated 



by lEkstrom et al. I (120081 ) who present evolutionary tracks for model 20 and GOM© stars. 
The tracks for the 2OM0 models are not immediately applicable since their portions that 
correspond to luminosity classes I and II occur in the cooler B-star domain (not included in 
our sample statistics). On the other hand, the tracks for 60Mq stars probably only apply 
to a small fraction of our sample of stars. However, the GOMq tracks do cover all luminos- 
ity classes, so we can use them to make a representative comparison of the predicted and 
obser ved changes in r otatio nal velocity with age. We plot in Figure 8 the 6OM0 models 



from lEkstrom et al. I (120081 ) that show the change in the average Vsini with age. Here we 
have multiplied the model equatorial velocities by a factor of n/A to represent the average 
projected velocity for random orientations. The thin solid lines show the velocity evolution 
for Galactic abundance stars for three values of the initial ratio of actual to critical angular 
velocity, and the thin dotted lines show the same for lower metallicity SMC stars. Above 
these are three tick marks indicating those ages associated with luminosities and temper- 
atures of luminosity class V, HI, and I stars according to the observational calibration of 



Martins et al.l (120051 ). For simplicity, we will assume that our sample groups of stars corre- 
spond to ages midway between 1 Myr (approximately when massive stars emerge from their 
natal, ultracompact H II region) and 1.83 Myr (the class HI age) for the unevolved group 
and midway between 1.83 Myr and the age at which they cool to T^g = 28 kK (the boundary 
between O- and B-supergiants) for the evolved group. We plot our derived average Vsini 
values for these groups (Table 6) versus their assigned mid-ages (relative to the GOM© tracks) 
as thick solid and dotted lines bounded by plus signs for the Galactic and SMC samples, 
respectively. We see that the modest declines observed are not too different from the model 
estimates for stars with relatively small rotational velocities. The tracks for lower mass stars 
are qualitatively similar, but their later, evolved portions where the rotational velocities 
change the most will be progressively shifted out of the 0-star domain into cooler classes, 
so that the predicted changes in mean rotational velocity between unevolved and evolved 
samples will be even smaller than shown in Figure 8. Thus, we emphasize that our simple 
comparison is probably only valid for considering relative trends with age and metallicity, 
and a detailed comparison would require additional tracks for other masses and sampling of 
the tracks based upon the observed mass and Vsini distributions. 

We also are interested in whether the initial V sin i distribution is the same in the three 
environments. In Figure 9, we plot the CDFs of all three dwarf samples. Statistically all three 
have p values that are above the cut-off of 5%, and certainly at the high V sini end all three 
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look extremely similar, indicating that the maximum rotational velocities are very similar. 
Below 200 km s~^, the behavior of the three distributions varies. Again we see that the shape 
of the CDF for the SMC dwarfs is dissimilar from that of the counterparts in the Galaxy and 
the LMC between 80 — 190 km s~^. Below 80 km s~^, the fraction of slowly rotating stars 
varies between the three environments, with the Galactic sample in between the SMC and 
LMC. In fact the largest divergence between the LMC and SMC CDFs comes a t V sin i ^ 65 



km s ^, resulting in a p value just above the statistically significant level. iHunter et al. 



(120081 ) showed an analogous plot, but for objects with M < 25Mq. Similarly they find good 
agreement between the LMC and Galactic samples, with the SMC sample lying slightly 
beneath the other two, especially in the region below 200 km s~^. It is interesting that there 
are no very slow rotators in our LMC dwarf sample. At these low velocities, the dominating 
effect must not depend upon metallicity since the Galaxy's metallicity, with significantly 
higher Z, is situated between that of the two low Z samples. We conclude that the initial 
velocity distribution in our three unevolved samples are statistically indistinguishable. 

Finally we compare the evolved stars in each environment to examine the relative frac- 
tions of stars with large Vsini (Fig. 10). Looking at the distributions at the high velocity 
end, we see that a slightly larger fraction of evolved stars in the SMC have Vsini values 
larger than 200 km s~^ than in the Galaxy or LMC, which might support theoretical predic- 
tions. But the large D values that are found between the LMC and SMC samples are again 
at the low velocity range, ~ 80 km s~^. The maximum differences between the Galaxy and 
both the LMC and SMC CDFs occur near this velocity, with both of the low metallicity 
environments having a larger fraction of stars with smaller values compared to the Milky 
Way sample. In fact the trend is very supportive of our earlier suggestion that metallicity 
plays an important role i n mac roturbulent broadening in evolved 0-type stars. A recent 



paper by ICantiello et al.l (120081 ) discusses the origin of atmospheric turbulence in massive 
stars by sub-surface convection zones that are driven by Fe-peak element ionizations. In 
their simulations, the threshold luminosity for the occurrence of an iron convective zone is 
ten times lower at Zmw than that for Zsmc- Our results support their prediction that tur- 
bulence will increase with metallicity. The SMC evolved sample has the largest fraction of 
slow rotators, followed by the LMC and then the Galaxy. We note that the D values for the 
SMC vs. Galaxy and LMC vs. Galaxy are extremely similar, but result in differing p values 
owing to the smaller SMC sample size. The large divergence at small V sin i values between 
the Galaxy and LMC, and possibly SMC, leads us to reject the hypothesis that they are 
both from the same parent distribution. Again we stress that the differences that cause this 
are not the fractions of stars with large V sini, but those with small values. 



In conclusion, we find some support for the new stellar interior model predictions that 
massive stars in lower metallicity environments will remain at almost constant rotation rates 
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throughout their MS hfetimes. But we also see that this effect, loss of angular momentum 
during the MS, is more subtle than previously reported even at the higher Galactic metallic- 
ity. We find that the ratio of mean Vsini in the evolved to unevolved samples is 0.84, 0.90, 
and 0.88 for the 0-stars of the Galaxy, LMC, and SMC, respectively. This is consistent with 
the general expectation of spin down and the specific expectation that the spin down rates 
are large in the Galaxy where mass loss rates are higher. Note that the ratio of 0.84 for 
the Galaxy is an upper limit, since the apparent projected velocities of the Galactic evolved 
stars are systematically higher because of larger macroturbulence. With the exception of the 
SMC unevolved sample, the largest divergence between sample CDFs occur at low Vsini 
values. We have suggested that metallicity may play an important role in the development 
of macroturbulence in the photospheres of the evolved, massive stars. Finally, we once again 
see that small sample sizes, particularly for the SMC, hamper our ability to establish the 
true distribution of equatorial rotational velocities. We urge the continued observational 
efforts to determine the Vsini values for a large number of massive stars in differing Z 
environments to address this pressing and important issue. 
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APPENDIX: NOTES ON INDIVIDUAL STARS 



AV 80 - This star is included in the CDF for SMC evolved stars in accord with the individual 



notes in 



Walborn et al.l (120001 ). who state that its absolute visual magnitude and the strength 



of its UV wind features correspond to a luminosity class of III. 



AV 177 - The CCF has a double-peaked appearance, however the velocity separation is not 
large enough for a double Gaussian fit. The V sini value listed is of both peaks together. 
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AV 3 32 - Known as a double-lined spectroscopic binary classified as WN3 + 06.5 I (IMassey et al 



20001 ) ■ Only one peak is seen in the CCFs, probably that of the 06.5 I component, which 



moves in Vr by 45 km s ^ between the observations taken 6 months apart. 



AV 476 - Classified as 02-3 V + companion by iMassey et al.l (120051 ) . we see two peaks, 
separated by K of 363 km s~^ 

NGC346-MPG324 - The peak in the two observations moves by 110 km s~^ and changes 
width. 

SK -65 21 - There are four observations taken in sequence over 15 hours. All CCFs show 
double peaks that appear to be near or at maximum velocity separation of 350 km s~^. The 
peaks first move slightly farther apart and then closer over the four observations. 



SK -66 172 - Classified as 02 III(f*) + OB (IWalborn et all 120021 ) . we do not see any sign of 
a companion in the CCF. 



SK -67 18 - iKoenigsberger et al.l (120031 ) classified this system as 03 -|- 08-BO I. We see both 
peaks in the CCF. 

SK -67 166 - Observed 63 times with FUSE over approximately a 5 week period, the peak 
does not move. 

SK -68 52 - A weak second peak is stronger in the CCF with AV327 than with SK -66 100. 
The peaks are separated by 350 km s^^. 

SK -69 94 - The FUSE spectrum o f S Dor is contam inated by the visual companion 13" 
away which is classified as a Bl I (jWolf et al.l Il980t ). The CCF here is almost certainly 
representative of this star, not S Dor. 



SK -69 223 - Classified WC4(+0B) bv iBartzakos et al.l (l200ll ). we do not see any sign of 
a second peak in the CCFs. However the main peak does move slightly between the two 
available observations taken 50 minutes apart. 



SK -70 91 - Classified a 02 III(f*) + OB bv IWalborn et al.l (j2002|), both stars have peaks 
in CCF. The peaks are separated by 255 km s~^. 

BI 173 - In the two observations taken two months apart, the main peak moves by 53 km 
s^^. There is a low level second peak, but it is difficult to distinguish it from the surrounding 
noise. 



HV 2241 - A known binary, we do not see any sign of the secondary, classified as early B 
jMokiem et al.ll2006h . in the CCF. 
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LMC X-4 - This spectrum is of the 08 III companion in this X-ray binary. 

PGMW 3209 - Only one extremely noisy observation is available from FUSE. We do not see 
any sign of a second peak in the CCF. 
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Table 1. Galactic Targets 
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Table 1 — Continued 
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58 


48 


Z9015101 


HD 115455 


07.5 IH((f)) 


69 


84 


85 


A1200701 


HD 116852 


09 HI 


136 


114 


95 


P1013801 


HD 122879 


BO la 


92 


69 


58 


B0710501 


HD 124314 


06 V(n){(f)) 


242 


285 


302 


P1026201 


HD 124979 


08 ((f)) 


220 


252 


280 


P1026301 


HD 148937 


06.5 f?p 


76 


77 


90 


C1680101 


HD 152233 


06 IH;(f)p 


104 


102 


102 


P1026702 


HD 152249 


OC9.5 lab 


99 


83 


72 


B0250401 


HD 152314 


09.5 HI-IV 


76 


86 


72 


P1026901 


HD 152405 


09.7 Ib-II 


77 


68 


56 


P3031801 


HD 152424 


OC9.7 la 


86 


74 


67 


Z9016301 


HD 152590 


07.5 V 


60 


76 


76 


B0710601-2 


HD 152623 


07 V(n)((f)) 


106 


97 


94 


P1027001 


HD 152723 


06.5 HI(f) 


111 


104 


107 


P1027102 


HD 154368 


09.5 lab 


102 


89 


78 


P1161901 


HD 157857 


06.5 HI(f) 


97 


99 


151 


P1027501 


HD 163758 


06.5 laf 


90 


84 


83 


P1015901 


HD 163892 


09 IV((n)) 


201 


287 


254 


P1027601-2 


HD 167659 


07 H(f) 


85 


77 


95 


P1028001 


HD 168076 


04 V((f)) 


98 


77 


116 


P1162201 


HD 168941 


09.5 II-IH 


116 


96 


82 


P1016501-2 


HD 175754 


08 n((f)) 


177 


207 


193 


P1016802 


HD 175876 


06.5 IH(n)(f) 


266 


274 


293 


P1016902 


HD 186890 


07.5 ni((f)) 


83 


88 


85 


P3030201 


HD 189957 


09.5 HI 


92 


92 


78 


P3030301 


HD 190429A 


04 If+ 


105 


77 


73 


P1028401 


HD 190864 


06.5 HI(f) 


88 


85 


91 


E0820501 


HD 191423 


09 IILri* 


436 


356 


336 


P3030101, E0821301 


HD 192639 


07 Ib(f) 


96 


110 


100 


C1710101, P1162401 


HD 193514 


07 Ib(f) 


94 


94 


94 


E0820701 


HD 199579 


06 V((f)) 


73 


70 


82 


P1162501 


HD 201345 


ON9 V 


91 


69 


58 


P1223001 


HD 207198 


09 Ib-II 


91 


71 


69 


P1162801 


HD 207538 


09.5 V 


51 


58 


36 


P1162902-3 


HD 210809 


09 lab 


117 


98 


81 


P1223101-3 


HD 216532 


08.5 V((n)) 


189 


225 


123 


A0510202 
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Table 1 — Continued 









< V sin i > 








Spectral 


Howarth et al. ( 


1997) AV 327 


SK -66 100 


FUSE 


Star 


Classification 


(km s-l) 


(km s^-"-) 


(km s-i) 


data set 


HD 217086 


07 Vn 


332 


310 


343 


E0820801 


HD 218915 


09.5 lab 


94 


75 


74 


P1018801 


HD 225094 


B3 la 


68 


77 




Z9120101 


HD 303308 


03 V({f)) 


111 


85 


137 


P1221601-2 


BD -48 3437 


Bl lab 


54 


44 




P1018401-3 


BD -60 2522 


06.5 (n){f)p 


178 


249 


228 


D1130101 


CPD -59 2600 


06 V((f)) 


142 


135 


182 


P1221401-3 


CPD -74 1569 


09.5 V 


64 


68 


59 


P1015301, U1092701 



Table 2. SMC Targets 



< V sin i > 





Spectral 




This paper 


Other 


Template 


FUSE 




Star 


Classification 


Rcf. 


(km s ) 


(km s ) 


Star 


data sets 


Notes 


AV 14 


03-4 V + Tioh 


I 


261 


212^ 


AV 327 


P1 1 75301 

-L J- -L 1 O'jyj J- 




AV 15 


Ofi 5 TTl'f") 


2 


95 


128^' 132^^ 135^ 

J_^0 _LO^ -LOvJ 


AV 327 


P1 1 501 01 




AV 26 


07 TTl-l-nph 


I 


141 


127'' 128'^ 


AV 327 


P1 1 76001 




AV 47 


08 iiicrfi'i 


2 


79 


75b j^ggc 


AV 327 


P1150202 




AV 61 


05 III 


3 


228 




SK -6fi ion 


F,5i 1 n2ni 




AV 69 




2 


94 


lOO'' 104*^ 


AV 327 


P1 1 5n3n3 




AV 70 


09 la 


3 


62 




AV 327 


A 1 1 80202-3 




AV 75 


05 ITI('f+~l 


2 


94 


109b iio<= 


AV 327 


P1150404 




AV 80 


04-fi nffin 

V_yrr U 111 1. I IJ 


2 


278 


324*= 270*^ 


AV 327 


E5111701 


* 


AV 83 


07 Iaf+ 


2 


71 


82*= 80^^ 70' 


AV 327 


P1176201 




AV 95 


07 TITCl'fl'l 


2 


76 


fiob QQc figa 

\Jk/ \J\D 


AV 327 


E054nini PiiRnsnR 




AV 96 


Bl la 


3 


70 




AV 327 


01 600101 


NO 


AV 104 


BO lb 


3 


56 




AV 327 


D1260101 




AV 119 


B2 II 


I 


48 




AV 327 


ZQ1 221 01 

ZJ X id ^ X U X 


NO 


AV 120 


09.5 III 


4 


70 




AV 327 


E5110301 




AV 175 


B2 lab 


3 


51 




AV 327 


Z9122001 


NO 


AV 177 


04 V 


3 


428 




SK -6fi ion 


nnnani 02 


*, NO 


AV 207 


07 V 


1 


75 




AV 327 


P 11 75901 

X XXI ^tJ\J±. 




AV 210 


B1.5 la 


3 


43 


37d 


AV 327 


D1620601 


NO 


AV 215 


BO la 


3 


92 


66^^ 


AV 327 


D1620401 




AV 216 


Bl III: 


3 


83 




AV 327 


D1620501 




AV 232 


07 Iaf+ 


5 


76 


106'' 74'' 


AV 327 


PI 030201 


SK 80 


AV 235 


BO law 


5 


90 




AV 327 


pinansni 


SK 82 


AV 238 


09.5 III 


5 


71 




AV 327 


P11766ni 

X XXI \J\J\J X 




AV 242 


BO 7 Taw 


5 


83 


54d 


AV 327 


Pn 76901 

X XXI \JtJ\J±. 


NO 


AV 255 


08 V 


g 


188 




AV 327 


E5110501 






Wo-0 1 


a 

D 


a I 




rt.V OZ ( 


HiOl iUDUl 




AV 264 


Bl I 


5 


116 




AV 327 


P1177001 


NO 


AV 266 


Bl I 


3 


61 




AV 327 


Z9122301 


NO 


AV 267 


08 V 


7 


268 




AV 327 


E5110701 




AV 321 


09 lb 


3 


357 




AV 327 


PI 150606 




AV 327 


09.7 Il-Ibw 


5 


69 


71*= 80^= 


AV 327 


P1176401 




AV 332 


06.5 I +WN3 


8 


90 




AV 327 


P1030401, X1050102 


*,SK 108, NO 


AV 372 


09 labw 


5 


98 


135'^ 


AV 327 


P1176501 




AV 374 


B2 lb 


9 


73 


47d 


AV 327 


C1600201 


NO 


AV 377 


06 V 


1 


51 




AV 327 


B1340101 




AV 378 


09 III 


3 


64 




AV 327 


P1150707 




AV 388 


04 V 


1 


147 


163'^ , 179'' 


AV 327 


P1175401 




AV 423 


09.5 II(n) 


5 


163 


186'' 


AV 327 


P1176701 




AV 435 


04 V 


6 


96 




AV 327 


E5111601 




AV 440 


07 V 


1 


62 




AV 327 


E5110901 




AV 446 


06.5 V 


1 


124 




AV 327 


E5111001 




AV 456 


09.5 lb 


4 


111 




AV 327 


Q1070101-4, 0106 


SK 143 


AV 462 


B1.5 la 


9 


57 


43'! 


AV 327 


Al 180301 


NO 
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Table 2 — Continued 



< V sin j > 





Spectral 




This paper 




Other 


Template 


FUSE 




Star 


Classification 


Ref. 


(km s"-*-) 




(km s"-') 


Star 


data sets 


Notes 


AV 469 


08 II 


1 


79 


SI'' 




AV 327 


P1176301 




AV 472 


B2 la 


9 


73 


33"^ 




AV 327 


C1600301 


NO 


AV 476 


02-3 V + comp 


10 


65(54) 






AV 327 


C0020301 


*, NO 


AV 488 


B0.5 law 


11 


87 






AV 327 


P1176801-3, P1030501 


SK 159 


NGC346-MPG324 


04 V((f)) 


5 


74,118 


70'= 


120'^ 105=" 70° 


AV 327 


P2030502, P1175601 


*, NO 


NGC346-MPG342 


05-6 V 


5 


58 






AV 327 


P2030401, P1175701 




NGC346-MPG355 


02 III(f*) 


5 


112 


112'= 120== 110== 


SK -66 100 


P1175201 




NGC346-MPG435 


04 III(n)(f) 


5 


93 






AV 327 


P2030201, P1175501 




SK 190 


08 laf 


13 


302 






AV 327 


E5111401 




SK 191 


B1.5 la 


9 


118 


92^ 




AV 327 


D1620201 


NO 


SK 197 


08.5 Ve 


12 


166 






AV 327 


E5111501 





Not e. — Spect ral Classificati ons from 
1997h , (4) lEvans et al. ( 2 0o3), (5) Iwa 



200C), (9) iDufton et al. 




iioi 



iGarmanv et alj (Il987h. (2 IWalborri et alj ll2000h. (31 ISrnith Neubig fc Bruhweiler 
born et alj bOO'j). f el lMassev et^ lll995l).(7') |CramDton fc GreaslevI h98^ . CSl lMassev et al. 



^ , ^ 

Massey et al. (2005), (11) 'Walborn ct al. (1995), (12) Court es et al.l lll995h. and ( 13) Massey 
2002D. Vsini valu es fr om (a)lMokiem et akl \2006). (b) Penny ct al. (2004), (c) Walborn ct al. (2000h. (d) iDufton et aLi I2006i1 . (e) 



iBouret et all ||2003|) . (fl lHillier et al.l ||2003| '). *Additional notes on individual stars are presented in the Appendix. 
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Table 3. LMC Targets 



< Vsini > 







Spectral 




This paper 


Other 


Template 


FUSE 






Star 


Classification 


Ref. 


(km s ) 


(km s ) 


Star 


data sets 


Notes 


SK 


-65 01 


BO. 5 I 


X 


74 




AV 327 


A04q0'^02 


NO 


SK 


-65 05 


09 II 


2 


153 




AV 327 


F,51 1 2S01 




SK 


-65 21 


OQ 7 Tab 


3 


57(58) 




AV 327 


P1020Q01-4 

L l.\J 1. ^ 


*, NO 


SK 


-65 22 


06 laf-l- 


3 


66 




AV 327 


PI 03 1002 




SK 


-65 44 


09 V 


4 


132 




AV 327 


P11733301-2 




SK 


-65 47 


04 If 


5 


208 




SK -fifi 100 

kJXV \J1J l.\J\J 


00020701 




SK 


-65 63 


09.7 I 


g 


77 




AV 327 


A0490701-2, B0861101, M1142001 




SK 


-66 01 


B1.5 la 


6 


78 




AV 327 


C1600601 


NO 


SK 


-66 18 




7 


100 


82* 


SK -66 100 


A0490101-2 




SK 


-66 78 


B1.5 I 


1 


65 




AV 327 


B1280301 


NO 


SK 


-66 97 


Bl lab 


8 


159 




SK -66 100 


E9570101 


Hen S35, NO 


SK 


-66 100 


06 11(f) 


3 


67 


84* 


AV 327 


PI 1 72303 




SK 


-66 106 


B2 la 


6 


78 




AV 327 


B1280901 


NO 


SK 


-66 118 


B2 la 


6 


167 




AV 327 


CI 600701 


NO 


SK 


-66 169 


OQ 7 la-l- 


3 


71 


72*^ 


AV 327 


PI 173801 


HD 269889 


SK 


-66 171 


09 la 


6 


90 




AV 327 


E5114102 




SK 


-66 172 




3 


101 




SK -fifi 100 

kJiV \J\J ±\J\J 


PI 1 72201 


*, NO 


SK 


-66 185 


BO lab 


2 


53 




AV 327 


A0490902 




SK 


-67 05 


09.7 lb 


3 


87 




AV 327 


PI 030703-4 




SK 


-67 14 


B1.5 la 


3 


88 




AV 327 


P1174201-3 


NO 


SK 


-67 18 


□ + O8-B0 I 


9 


80(67) 




SK -66 100 


nnqsosoi 


*. NO 


SK 


-67 28 


BO. 7 la 


Q 


76 




AV 327 


A0490201-2 


NO 


SK 


-67 38 


08.5 V 


7 


100 




AV 327 


E5112601 




SK 


-67 46 


B1.5 I 


1 


71 




AV 327 


A0491501 


NO 


SK 


-67 69 


04 ni{f) 


3 


121 




SK -fifi 100 


PI 171 701 1703 

X 1. ± 1 ± 1 \J 1. ^ X( XJtJ 




SK 


-67 76 


BO la 


11 


80 




AV 327 


P1031201 




OJA. 


-D ^ iuu 


oi ia 


a 
D 








■Rl OQflzini 


FTFl 9fiQPinzl l\TO 


SK 


-67 101 


08 II((f)) 


3 


103 


lOl'' 


AV 327 


P1173401-3 




SK 


-67 106 


BO I 


12 


99 


135'= 


AV 327 


AlllOlOl 




SK 


-67 107 


BO I 


13 


82 


103'' 


AV 327 


A1110201 




SK 


-67 111 


06 Ia(n)fp vax 


3 


259 


209'' 


AV 327 


C1550208-9, C1550214-5 




SK 


-67 118 


07 V 


7 


61 




AV 327 


E51 13401 




SK 


-67 119 


07 Ill(f) 


7 


235 




AV 327 


E5113501 




SK 


-67 150 


Bl la 


4 


66 




AV 327 


B1281001 


NO 


SK 


-67 166 


04 If-f 


3 


97 


157'' , 97* 


SK -66 100 


A1330101-162, A1330180 


* 


SK 


-67 167 


04 Inf-l- 


3 


155 




SK -66 100 


P1171901-2 




SK 


-67 169 


Bl la 


3 


62 




AV 327 


P1031603 


NO 


SK 


-67 174 


08 V 


7 


72 




AV 327 


D0980101 




SK 


-67 176 


07 Ib(f) 


6 


54 




AV 327 


D0980201 




SK 


-67 181 


B0.5 I 


11 


99 




AV 327 


B0900101 


NO 


SK 


-67 191 


08 V 


2 


100 




AV 327 


P1173101-2 




SK 


-67 256 


Bl la 


6 


66 




AV 327 


B1280201 


NO 


SK 


-67 266 


Ofpe/WN9 


14 


113 


144b 


SK -66 100 


B0900201 




SK 


-68 03 


09 I 


2 


70 




AV 327 


A0490401-2 
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Table 3 — Continued 



< V sin i > 





Spectral 




This paper 


WTiner 


Template 








Classification 


rtei. 


t^KlIl S J 


I^Km s ) 


otar 


data sets 


w otes 


SK -68 16 


07 III 


2 


100 




AV 327 


E5112201 




SK -68 41 


BO. 5 la 


6 


77 




AV 327 


P1174101-2 


NO 


SK -68 52 


BO la 


3 


70(55) 




AV 327 


P1174001 


*, NO 


SK -68 75 


BI I 


1 


78 




AV 327 


A0490501-2 


NO 


SK -68 111 


BNl la 


15 


72 




AV 327 


B1280501 


HD 269668, NO 


SK -68 129 


BI I 


5 


96 




AV 327 


B0860501 


NO 


SK -68 135 


ON9.7 Ia+ 


16 


68 


70'' 


AV 327 


P1173901 


HD 269896 


SK -68 137 


03 III(f*) 


16 


107 




SK -66 100 


E5114201 




SK -68 155 


BO. 5 


11 


96 




SK -66 100 


B0860701 


NO 


SK -68 171 


BI la 


6 


71 




AV 327 


A0490801 


NO 


SK -69 43 


BO. 5 la 


6 


102 




AV 327 


B1280601 


HD 268809, NO 


SK -69 50 


07 If 


7 


178 




AV 327 


E5112001 




SK -69 52 


B2 la 


6 


105 




AV 327 


P1174301 


NO 


SK -69 59 


BO la 


11 


111 




AV 327 


P1031103 




SK -69 94 


AOc/LBV 




127 




AV 327 


Z9050201 


*,S Dor = HD 35343, NO 


SK -69 104 


06 Ib(f) 


3 


117 


97'' 


AV 327 


P1172401 


HD 269357 


SK -69 124 


09 lb 


3 


96 




AV 327 


P1173601-2 




SK -69 223 


WC4(-fOB) 


17 


84 




AV 327 


P2150201-2 


*HD 38029, Brey 67, NO 


SK -69 237 


BI la 


15 


86 




AV 327 


B1280701 


NO 


SK -69 279 


09 f 


2 


84 




AV 327 


B0860801 




SK -70 13 


09 V 


7 


84 




AV 327 


D0981001 




SK -70 32 


09.5 II 


18 


97 




AV 327 


E5112501 




SK -70 60 


04-5 V:n 


13 


256 




SK -66 100 


P1172001 




SK -70 69 


05 V 


3 


98 


131=- 


AV 327 


P1172101 




SK -70 78 


B0.7 la 


6 


107 




AV 327 


B1280801 


HD 269074 


SK -70 79 


BO III 


7 


86 




AV 327 


B0900301 




SK -70 85 


BO I 


1 


63 




AV 327 


A0491301, A0491303 




SK -70 91 


02 III{f*)-)-OB 


3 


99(78) 




SK -66 100 


P1172501 


*, NO 


SK -70 115 


07 lb 


4 


164 




AV 327 


P1172601 




SK -70 120 


BI la 


6 


75 




AV 327 


A0491001-2 


NO 


SK -71 08 


09 II 


2 


75 




AV 327 


A0491401-2 




SK -71 19 


07 V 


4 


242 




AV 327 


E5113001 




SK -71 41 


08.5 I 


7 


53 




AV 327 


D0980501 




SK -71 45 


04-5 Ill(f) 


3 


117 


125'' 


SK -66 100 


P1031501-4 


HD 269676 


SK -71 46 


04 If 


7 


79 




SK -66 100 


D0980601 




SK -71 50 


06.5 III 


2 


296 




AV 327 


A0491201 




BI 12 


06 III 


4 


116 




AV 327 


E5111802 




BI 13 


06.5 V 


2 


118 




SK -66 100 


E5111901 




BI 130 


08.5 V((f)) 


7 


95 




AV 327 


E51 12902 




BI 170 


09.5 lb 


3 


84 




AV 327 


P1173701 




BI 173 


08 11: 


3 


116 




SK -66 100 


P1173201-2 


*, NO 


BI 208 


07 V 


2 


196 




AV 327 


P1172702-5 




BI 237 


02 V((f*)) 


10 


83 


126* 


SK -66 100 


E5114001 




BI 253 


02 V((f*)) 


10 


155 


19P 


SK -66 100 


C0020601-2 
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Table 3 — Continued 



< V sin j > 





Spectral 




This paper 


Other 


Template 


FUSE 




Star 


Classification 


Rcf. 


(km s"-') 


(km s^-^) 


Star 


data sets 


Notes 


BI 272 


07 III 


19 


212 




AV 327 


P1172901-2 




HV 2241 


07 III + BO III 


20 


71 




AV 327 


F9270603 


*, NO 


LH47-338 


07 V((f)) 


5 


99 




AV 327 


P2320101 


N44C Star 2 


LMC X-4 


08 III 


21 


322 




AV 327 


A0440102 


*, NO 


LMC2-703 


09.5 V 


7 


95 




AV 327 


D0981401 


D301-1005 


LMC2-755 


08 V 


7 


57 




AV 327 


D0981501 


D301-NW8 


PGMW 1486 


06.5 V 


22 


93 




SK -66 100 


D0300101 




PGMW 3053 


05.5 I-III(f) 


22 


110 




SK -66 100 


BOIOOIOI 




PGMW 3070 


06 V 


22 


72 




AV 327 


B0100201 




PGMW 3073 


06.5 V 


22 


71 




AV 327 


D0300201 




PGMW 3102 


07 V 


22 


68 




AV 327 


D0300302 




PGMW 3120 


05.5 V{(f*)) 


22 


111 


95*^ 


SK -66 100 


B0100301 




PGMW 3157 


BCl la 


22 


88 




AV 327 


B0100401 


NO 


PGMW 3168 


07 11(f) 


22 


78 


96a 


SK -66 100 


B0100501 




PGMW 3204 


06-7 V 


22 


78 


130^ 


SK -66 100 


D0300402 




PGMW 3209 


03 III(f*)+OB 


22 


166 




AV 327 


B0100601 


*, NO 


PGMW 3223 


08.5 IV 


22 


147 




SK -66 100 


B0100701 




PGMW 3264 


03-6 V 


22 


56 




SK -66 100 


B0100801 





Note. — Spectral Classifications fro m ( I'lljaxon et all ll200j) , f2 )IConti et all lll986h . ( 3) IWalborn et al.l ll2002f) . 
f4llSmi th Neubig fc Bru hweiled lll999lv (5')lM'assev et al.l ||2000|'I . (6)[Fitzpatrickl (Il98^. (7VMassev et al. | lll995h7(8 ) 
Gummc rsbach et al. (199^ . (9) |Koenigsberger et al.l2003l), (lO'llMasscv et al. (2005), (11) Rousseau c t alj lll978ll 
(12) [pihren bach & Duflot iwsj ) , (13) iRobert et^ ||2003^ , (14) Boliannan & Walborn (19^ 
1991 ). (16) ,Walborn (19771), (1 7 )lBartzakos et alJ fcOOll). (18)|Massev (2002). (19)£_lantc (199 



Walborn (1989), f 15 ) iFitzpatrick 



20011) , (21) Negueruela fc^Coel (l2002h . and (22) IParker et al.l ir992h . Vsini values from (a) 



and (b jPennv et al.l (l2004h ~ *Additional notes on individual stars are presented in the Appendix 



(20)IOstrov et al. 
Mokiem et al.l (|2007^ 
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Table 4. SMC Stars with Previous ^ sini Values 





Spectral 


< Vsinj > 




Star 


Classification 


(km s~-^) 


Reference 


346-010 


07 IIIn((f)) 


313 


M06 


346-012 


Bl lb 


29 


M06 


346-018 


09.5 Ille 


138 


M06 


346-022 


09 V 


55 


M06 


346-025 


09 V 


138 


M06 


346-031 


08 Vz 


18 


M06 


346-033 


08 V 


188 


M06 


346-046 


07 Vn 


340 


M06 


346-050 


08 Vn 


357 


M06 


346-051 


07 Vz 


18 


M06 


346-066 


09.5 V 


129 


M06 


346-077 


09 V 


177 


M06 


346-090 


09.5 V 


188 


M06 


346-093 


BO V 


187 


M06 


346-097 


09 V 


22 


M06 


346-107 


09.5 V 


55 


M06 


346-112 


09.5 V 


143 


M06 


330-013 


08.5 III((f)) 


73 


M06 


330-052 


08.5 Vn 


291 


M06 


NGC346-MPG368 


04-5 V({f)) 


55 


P04 


NGC346-MPG113 


OC6 Vz 


<40 


P04 


AV 243 


06 V 


62 


P04 


NGC346-MPG487 


06.5 V 


65 


P04 


AV 220 


06.5f?p 


<40 


P04 


NGC346-MPG682 


08 V 


71 


P04 


NGC346-MPG12 


O9.5-B0 V 


67 


P04 


AV 170 


09.7 III 


54 


P04 



Not e. — M06 = iMokiem et~al] ll2006h and P04 = IPenny et al.l 
i2004h 
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Table 5. LMC Stars with Previous Vsini Values 





opcctrcil 






star 


Cl&ssificSytrion 


(km s ) 


Jrx.GlGrGIlCS 


Nll-004 


OC9.7 lb 


81 


M07 


Nll-008 


B0.7 la 


83 


M07 


Nll-026 


02 III(f*) 


109 


M07 


Nll-029 


09.7 lb 


77 


M07 


Nll-031 


ON2 III(f*) 


116 


M07 


Nll-033 


BO Illn 


256 


M07 


N 11-036 


B0.5 lb 


54 


M07 


N 11-038 


05 II{f-f ) 


145 


M07 


Nll-042 


BO III 


42 


M07 


Nll-045 


09 III 


105 


M07 


Nll-051 


05 Vn((f)) 


333 


M07 


Nll-058 


05.5 V((f)) 


85 


M07 


Nll-060 


03 V{(f*)) 


106 


M07 


Nll-061 


09 V 


87 


M07 


Nll-065 


06.5 V((f)) 


83 


M07 


Nll-066 


07 V((f)) 


71 


M07 


Nll-068 


07 V((f)) 


54 


M07 


Nll-072 


B0.2 III 


14 


M07 


N 11-087 


09.5 Vn 


276 


M07 


Nll-123 


09.5 V 


110 


M07 


HD 269810 


02 III(f*) 


173 


P04 


SK -69 212 


05 Ill(f) 


210 


P04 


SK -67 51 


06.5 III 


105 


P04 


CD -68 264 


08 V 


139 


P04 


LH58-52A 


08 II 


110 


P04 


LH58-19A 


09 II 


178 


P04 



Note. — M07 = iMokiem et"al] l|2007h and P04 = 
IPennv et all JSooi) 
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Table 6. Kolmogorov-Smirnov Statistics 





"1 


< V sin i > 


"2 


< V sin i > 




p 


Sample 1 


stars 


(km s ^ 


) Sample 2 


stars 


(km s-l) 


D 


(%) 




Uncvolvcd Vs. 


Evolved Within Each Environment 






Galactic uncvolvcd 


79 


129.8 


Galactic evolved 


56 


109.5 


0.296 


0.5 


LMC unevolved 


36 


132.0 


LMC evolved 


42 


118.6 


0.151 


74 


SMC unevolved 


32 


116.2 


SMC evolved 


19 


101.8 


0.289 


23 








Unevolved Samples 










LMC unevolved 


36 


132.0 


Galactic unevolved 


79 


129.8 


0.180 


37 


SMC unevolved 


32 


116.2 


Galactic unevolved 


79 


129.8 


0.172 


47 


SMC unevolved 


32 


116.2 


LMC unevolved 


36 


132.0 


0.295 


8 








Evolved Samples 










LMC evolved 


42 


118.6 


Galactic evolved 


56 


109.5 


0.321 


1 


SMC evolved 


19 


101.8 


Galactic evolved 


56 


109.5 


0.297 


13 


SAIC evolved 


19 


11)1.8 


L.MC evolved 


12 


IIS.O 


i).i:-!7 


95 
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Fig. 1. — Evolutionary tracks from ISchaller et al.l (119921 ) with an overlai n lines at Tpff 
26000 K and T^g = 14900K. The first corresponds to spectral class 09.7 I by lCrowther et al. 



( I2OO2I ). presumably the lower limit of effective temperatures to be assigned to O spectral 
type stars. All stars to the right of this line will be typed as B types or cool er. The second 
overl aid line corresponds to the temperature for spectral classification B3 I (iBohm-Vitense 



198ll ). Stars with this temperature may have progenitors with initial masses as low as 7Mq. 
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Fig. 2. — The cross-correlation functions and Gaussian fits for the Galactic stars (a) 
HD 46202 (09 V; Vsini = 37 km s"^), b) HD 152233 (06 III:(f)p; Vsmi = 104 km s"^), (c) 
HD 96917 (08.5 Ib(f); Vsmi= 176 km s"^), (d) HD 41161 (08 Vn; Vsini ^ 304 km s"^). 
The bumpy nature of the ccf for HD 41161 may be indicative of non-radial pulsations. 




300 



Fig. 3. — The projected rotational velocities of iPennyl (119961 ) or iHowarth et al.l (119971 ) from 
lUE spectra vs. the Gaussian width of the CCF's from the template (a) AV 327 (09.5 
Il-Ibw) and (b) SK -66 100 (06 11(f)). 




V sin i (Mokiem et al. 2006,2007)(km s"^) 




50 100 150 200 250 300 350 
V sin i (Penny et al. 2004)(km s"') 



Fig. 4. — Proj ected rotatio n al vel ocities from Tables 2 and 3 versus (a) iMokiem et al.l (12006 



20071 ) and (b) iPenny et al.l (120041 ) for those targets in common. 
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Fig. 5. — Cumulative distribution functions of V sini values for unevolved and evolved stars 
in the Milky Way. 
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Fig. 6. — Cumulative distribution functions of V sini values for unevolved and evolved stars 
in the LMC. 
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Fig. 7. — Cumulative distribution functions of V sin i values for unevolved and evolved stars 
in the SMC. 
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Fig. 8. — Predicted changes in mean V sini with age for 60Mq models from ?. The thin 
sohd hnes show the variations for Galactic abundance stars with initial ratios of equatorial 
to critical angular velocity of 0.5, 0.3, and 0.1 {from top to bottom), while the thin dotted 
hnes show the same for SMC abundance stars. All these tracks are terminated at an age 
where the star cools enough to become a B-supergiant. The thick sohd and dotted lines that 
connect plus signs show the observed change in average Vsini between the unevolved and 
evolved groups for the Galactic and SMC samples, respectively. 
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Fig. 9. — Cumulative distribution functions of V sini values for unevolved (luminosity classes 
IV - V) stars in the SMC, LMC and Galaxy. 



-40- 




Fig. 10. — Same as Figure 9 but for evolved (luminosity classes 1 -11) stars in the same three 
environments. 



